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. Abstract. The sound propagation is usually described by a linear homogeneous wave equation, 

' though the air flow in a duct is described by the gas dynamics equations, using a variable cross 

section, which corresponds to a non linear non homogneous system. The aim of this paper is to 
exhibit a common periodic solution to both models, with several free parameters such as frequency 
, ^ . or amplitude, able to represent any sound. By taking in account a friction term linked to the material 

(wood or brass for instance) of the duct, it is possible to build an analytic such solution when the 
cross section fuUfills some condition which corresponds exactly to the general shape of the wind 
instruments. The conclusion is that in the wind intrumcnts, the shape yields linearity. 
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P^'. 1. The gas dynamics models with source terms 

I The duct is supposed to have a cyhndrical symetry along its axis which can be a straight hne (as 
• for a clarinet) or a curve line (as for a horn) and the abscissa along this axis is denoted by x > 0, 
. with X = corresponding to the narrowest side of the duct. The variable cross section is denoted 
by a{x) > 0, the air density and velocity are p and u, both depending on the position x and time 
t. Since people blows at x = 0, the velocity u is expected to be positive. The mass conservation is 
described by 

a{x) pt + (a(x) p u)^ = , 
' which becomes the transport equation 

cn ; (1.1) + = , 

by using the transported quantity q = a{x)p and the flux m = a{x)pu. This flux is ruled by the 
0\ ', general dynamic equation 

(1.2) mt + 2 u + (c^ - u^) + S{q,m) = 0, 

I where c = c(g) is the sound velocity and S" is a source term to be detailed later. The sys- 
^ ■ tem (11.11) . (11. 2p is an inhomogeneous form of the isentropic Euler equations, with the usual wave 
velocities u — c and u + c. 

A traveling-wave is expected to have the form q{x, t) = q{x — At), m = m{x — At) , with some 
constant A > corresponding to a velocity, so that both q and m are solution to the single wave 
equation: qt + A q^ = , mt + A = . 

Following the notion of source waves, such that Roll waves in hydraulics (see [1] or [3]), such 
solutions are exhibited by writing m = m{q) in (11. ip and (II. 2p . One gets 

({m'{q) -u)^ -c^)q^ = S{q,m{q)) . 
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Next, by introducing a function ip{q) such that 

{■m'{q) — — (? 



ip'{q) 



S{q,m{q)) 
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we obtain the two equations tp'{q) Qx = 1 and 'ip'il) Qt = —m'{q) ip'{q) Qx = —m'{q) . 
By integrating the first equation with respect to x we get 'ip{q) = x — K{t) , where K{t) does not 
depend on x. Next, a derivation with respect to t , using 'ip\q) qt = —m\q) , leads to m\q) = K'{t) . 
Another derivation, with respect to x this time gives m"{q) qx = . Since the expected wave is not 
flat, we select m'\q) = 0, which leads to the two expressions m{q) = A q — B , K'(t) = A , 
with some constants A and B , and the relation 

(1.3) -^{q) = x-At, 

with the integration constant put inside the expression of ip{q). From (11. 3p it is clear that q and 
then m = Aq — B are both functions of x — At , and therefore solutions of a linear wave equation 
with the constant velocity A (traveling- waves) . The expression of ip'{q) gets simpler by noticing 
that u = — = A — — so that m'iq) — u = — , which leads to 

q q q ' 

R2 _ „2 2 

(1.4) i,'{q) 



q' Siq,Aq-B) 

It is important to notice that the hypothesis S{q, m) ^ is fondamental for this result, which does 
not apply for homogeneous systems. Another fondamental hypothesis is that ip'{q) is a function of 
q only. The result does not apply to the case of a source term depending on q, m and x for instance 
or to a sound speed depending on q and x. In the next sections we look for state laws and duct 
shapes for which these hypotheses are verified, and conclude that it is genuinely the case for the 
usual shapes of wind instruments, as fiutes, clarinets or horns. 

2. Application to an air flow in a duct 

The cross section a(x) is supposed to be increasing with x (a'(x) > 0). We look for expressions of 
c(g) and 5'(g, m) compatible with the previous hypotheses, that is independant of x. The expression 
of c(g) comes from the pressure P using = ^ , where 7 = 1.4 is the adiabatic constant, and 
the pressure is linked to the density p and the temperature T (in ° Kelvin) by a state law. A well 
known state law is given by the Boyle-Mariotte law P = KqpT, with Kq = 287.06 in M.K.S. units, 
with P expressed in Pascals, which provides the following sound velocity 

c = ^J^^{K^ . 

Since the usual temperatures are greater than 270 " K, the variations of the profile of c are negligible, 
and we can fix the temperature to a reference value Tq and take c = cq = Kq Tq . This 
corresponds to the " isothermaT hypothesis. 

Another state law is given by the " isentropic" case P = K p'^, with K = 69259.5 in M.K.S. 

7-1 

units. The expression of the sound speed reads c = ^"^K , which is not a function of q only. 
Since (? = P'{p), we can derive the expression of the source term at rest, that is when the velocity 
is zero and the pressure is constant: u = 0, ^ = . The dynamical equation (11.21) reduces to 
c^qx + S{q, 0) = . Since qx = a{x) px + a'{x) p , we obtain a'{x) p + S{q, 0) = , which 
leads to two possible of the source term: 



S{q,0) 



— Cq ^^y^^ q in the isothermal case 



-7 K q"^ in the isentropic case . 



^0 a{x) 
a' 

alx)-' 

We introduce a friction term, of the Strickler type, to select the following source term 

S{q,m) = S{q,0) + k \u\u , 
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where k is the Strickler friction (dimensionless) coefficient. 

The retained model is a combination between the isothermal and the isentropic cases. The 
sound velocity is a constant c = cq = ^/'jKqTq, with a fixed temperature Tq. For instance with 
To = 300i^ we get Cq = 347.225 ms-\ 

The selected source term is S{q,m) = k\u\u — jKq'^ ,which becomes independent on x if 
we impose a profile such that 



a(x) L)^ 
(2.1) — = Constant = k — , 

a{x)' jK 

where D is a constant. Expression (12.11) is a differential equation whose solutions are 



alx) 



(7-1) k (xo-x) 

where Xq is a constant to be taken larger than L the length of the wind intrument. The radius of 



the duct is given by r{x) = y , and we recognize the usual shape of the wind instrument as 

seen on the figure below, where kD^ = 75 'jK . The value of Xq corresponds to an infinite radius, 
or to a maximal length for a wind instrument. By the way, using the isothermal case leads to 
another expression of the cross section which reads a{x) = exp (Cx) where C and Gq are positive 
constants, and allows an infinite length. 




Figure 1. The shape and the phase plane 



The dynamic equation (11.21) becomes 
(2.2) mt + 2um^ + (cq - u^) + k {u^ - D^q^) = , 

where the source term is independent on x which allows waves of the form m = 



Aq-B. 
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3. The wave profile 

The updated form of ip^q) reads 
(3.1) ^'(g) - 



k (^{Aq - Bf - q 



2 ^7 + 2 



and we consider a so-called Reference state M* = (g*, m,,,) , with > 0, m* > and set m* = ^ and 
take A = u^: + Cq and B = cq. The value g,,, is a possible realistic value ofg which is a root of 
the numerator in (13. ip . Since 4''{q)qx = 1, either g^; becomes infinite when g = g,,, or g* is also 

a root of the denaminator in (13. ip and in this case = D g* , and belongs to the set 
5*0 = {(g, m) I S{q, m) = } in the phase plane (g, m). The line m = (u^, + cq) g — g*Co cuts 5*0 in 
M*, coming from the set {S' < 0} for g < g^, and going towards {S" > 0} for g > g*. Since 

^'(g) - - ^° 



g2 5 (g, (m^ + Co) g - g*co) 

we see that for ilj'{q) < for g < g*, 5 < and for q > q*, S > 0. Since ip'{q)qx = 1, we get 
that q^ < 0, which means that the wave profile corresponding to the line m = {u^ + Cq) g — g*Co is 
decreasing. 

Now let us consider in {S* > 0} a point Mi = (gi, mi) on the line m = (u* + cq) g — g*Co ,and 
then the line whose reference point is Mi, that is the line of equation m = {ui + cq) g — giCo. Along 
this line, we have 

^'(,) = ^L^^!l^g ^ 

g2 5" (g, (iti + Co) g - giCo) 

which is positive for g < gi. This means that the corresponding profile is increasing. By the same 
way, let us consider in {S < 0} a point M2 = (g2, ^2) on the line m = (m* + cq) g — g^co ,and then 
the line whose reference point is M2, of equation m = {u2 + Cq) q — g2Co. Along this line, we have 



g2 S (g, {u2 + Co) g - g2Co) 

which is positive for q > q2- This means that the corresponding profile is increasing too. This is 
reprensented in a window on both figure, above and below. Up to now, the two points Mi and M2 
have been chosen without any restriction, and it is possible to get the two lines m = {ui + cq) g — giCo 
and m = {u2 + cq) g — g2Co cutting 5*0 at the same point Mq G 5*0. This can be done by choicing 
first Mo e 5*0 with go > g*, then searching a point M = (g, {u^ + Co)g — g*Co) such that Mq belongs 
to the lines of slope Uj — cq passing through Mj, j = 1,2. This corresponds to the equation 

(3.2) g+M^l = L + 2co - Z^go2 ) ^ , 

g \ I Co 

which have two solutions which are actually gi and g2 with g2 < g* < go < gi as shown in the figure 
below. 



The profiles are computed by integrating ip'il) g^ = 1 for each part of the profile, using a 
numerical integration method, with 
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Figure 2. The profile of the wave 



/// \ _ ^0 

^^'^*' ~ k{2{u^ + c,)u,ql-D^{^ + 2)q:^^) 

for q = q^. We first compute Mi and M2 from fl3.2p and starting from Mi we compute the two 
branches MiMq and MiM^,, next starting from M2 we compute the two branches M2M0 and M2M*. 

4. The wave propagation 

The profile is made of a sequence of 3 travelling waves: a rear part from the state Mq to Mi with 
the constant velocity Ai = ui + cq = + 2co — co|^ , a central part from the state Mi to M^ and 
M2 with the constant velocity A^, = u^, + cq and a front part from the state M2 back to Mq with the 
velocity A2 = M2 + Cq = m* + 2co ~ . Since q2 < q* < qo < qi, we have Ai > A^: > A2 ■ When 
and go are near, which is the case in reality, these velocities are also very near and we can consider 
that the whole wave travels at the velocity + cq. The figure above corresponds to a choice of the 
parameters such that g* and go are not too near and then the three lines in the phase plane do not 
seem superimposed for the reader. 

The effect of the different velocities together with the mass conservation principle is the emer- 
gence of two small shock waves, one on the top of the rear wave, in Mi and another at the lowest 
part of the front wave in M2. The amplitude of these two shocks increases during the propagation 
but stay imperceptible since the travel time in the duct is very short, of the order of 10~^ s. How- 
ever, such mass effects can be observed in the real world, with similar equations, sometimes at a 
very different scale in space and time (see for instance [4], on the modelling of the Rogue Waves 
using a small difference of the velocities of two long waves on the ocean). 

The sound production is produced by the concatenation of such waves with the same parameters, 
involving a constant wavelength and therefore a constant frequency. 

5. Some remarks and conclusion 

The wave profile presents a decreasing part, from Mi to M2, which never degenerates into a shock 
as often expected in hydrodynamics (as for the Roll waves, see [1] or [3]). This comes obviously 
from the source term. However, the occurrence of a shock may be analysed. Since the sound speed 
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is a constant, the pressure law corresponds to P{q) — c^q and the shock condition derived from the 
well known Rankine-Hugoniot relations for a shock wave linking Mi to M2 read 



U2-U1 / P{q2) - P{qi) Co 



q2-q\ V ^2 - ?2 V?i?2 

We shall compute the rate by two other ways. First we write that Mq belongs to the hues 

issued from Mi and M2, with the respective slopes Ui + Cq and 1*2 + Cq, that is 

(mo = ) (mi + Cq) go - coqi = {u2 + Co) go - cog2 , 

which gives 

U2 - Ui _ Co 

q2 -qi qo ' 
Next we write that M*, Mi ans M2 are aligned or 

lil = li* + Co — Co — , U2 — U* + Co — Co — , 

qi Q2 

U2 - Ui coq* 



which gives 



?2 - qi gi?2 



We have obtained 



U2-U1 \ 1 1 q* ■, , — 

Co = Co — = Co , hence q^ = V?i?2 = ?o ■ 



92 - gi / ^0 gig2 

Thus = Uq and {u^ =)ui + Cq — Cq^ = Mi — Co + Co^, which gives 2g*gi — qf — ql = 0, that is 

(g* — gi)^ = 0, or qi = q*. The same arguments give q2 = q*- Our hypothetical shock is reduced to 
one point, which means obviously that no shock occurs. 

The wave velocity A is close to Co, but a priori different since A — cq — U2 for the front part of 
the wave for instance. Depending on the (still realistic) values of the parameters of the model, U2 
can be negative or positive, and a possible constraint on the parameters may consist into setting 
U2 = 0. The arguments for this approach are settled in the reference [2]. 

This construction of a wave should also feed the discussion between linear or non linear models 
(see for instance [5]). It seems, as a conclusion, that the difference is not so important when a source 
term is present. By taking in account a non linear diffusion term instead of a friction term, other 
travelling waves may be expected. The main conclusion is that the linearity is naturally induced 
by the shape in a wind instrument. 

This paper is also available on the website "conservation laws preprint server", 2008-010. 
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